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We investigate how the different velocities characterizing the low-energy spectral properties and 
the low-temperature thermodynamics of one-dimensional correlated electron systems (Luttinger 
liquids) affect the transport properties of ring-like conductors. The Luttinger liquid ring is coupled 
to two noninteracting leads and pierced by a magnetic flux. We study the flux dependence of 
the linear conductance. It shows a dip structure which is governed by the interaction dependent 
velocities. Our work extends an earlier study which was restricted to rather specific choices of the 
interaction parameters. We show that for generic repulsive two-particle interactions the number of 
dips can be estimated from the ratio of the charge current velocity and the spin velocity. In addition, 
we clarify the range of validity of the central approximation underlying the earlier study. 

PACS numbers: 71.10.Pm, 72.10.-d, 73.63.Nm 
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I. INTRODUCTION 

In the presence of a two-particle interaction the low- 
temperature thermodynamics and the low-energy spec- 
tral properties of a wide class of one-dimensional (Id) 
electron systems is described by the Luttinger liquid (LL) 
phenomenology^ One of the characterizing properties of 
a LL is the complete separation of the fundamental col- 
lective spin and charge excitations. In the low-energy 
limit and for a fixed number of left and right moving 
electrons the Hamiltonian can be written as the sum of 
two terms describing free bosons for the spin and charge 
degrees of freedom, both having a linear dispersion with 
velocities v s and v c . Furthermore, in Id the bosonic de- 
grees of freedom span the entire low-energy Hilbert space 
at fixed particle number. This has to be contrasted to 
the collective spin and charge excitations of higher di- 
mensional Fermi liquids which only give a certain part 
of the spectrum. In many theoretical studies the effect 
of spin-charge separation was discussedi^^ and there 
have been several attempts to find experimental indica- 
tions for it£ 

The low-energy physics of a LL is parameterized by 
four (in general) independent and interaction dependent 
fundamental velocities vn, c , Vn,s, vj iC , and vj <s . The 
first two are relevant when the total charge or spin are 
changed, while the later two are the velocities character- 
izing the charge and spin current. From these, v s and v c 
can be computed as^ 



c/s — y/VN^c/s Vj,c/s 



(1) 



Here we are interested in the role of spin-charge sep- 
aration or more generally in the role of the different in- 
teraction dependent velocities on transport through LLs. 
There is no specific effect on the linear conductance G of 
a linear LL wire connected to a Fermi liquid source and 
drain. A different situation occurs if a closed LL ring 
coupled via tunnel barriers to two noninteracting leads 



is considered. This geometry allows for interference and 
an additional parameter can be added to the problem, 
namely a magnetic flux <j) piercing the ring. The linear 
conductance then becomes a periodic function of <f> with a 
periodicity of one flux quantum cf>Q. We here chose units 
such that 0o = 27r. 

In Ref. it was argued that for v c /v s — q/p, with q and 
p being "small" odd numbers, G((f>) shows characteristic 
dips whose number for <j> £ [0, 2ir] is given by qvp/v c = 
pvp/v s , with vf being the noninteracting Fermi velocity. 
The appearance of such dips might thus be interpreted 
as an indication of spin-charge separation. Restricting 
the ratio of the charge and spin velocities to such cases, 
however requires rather specific interaction parameters 
which will be hard to realize in experimental systems. 
The authors study a continuum model, the Tomonaga- 
Luttinger model (TLM). The central approximation of 
this work — which can be applied for the above given ratio 
of v c and v s — is related to the idea that an electron can 
pass the ring only if the charge and spin can "recombinc" 
at the right (and left) contact (for details see below). A 
dip structure of G{<j>) was also obtained for a small ring 
described by a lattice model with strong local Coulomb 
interaction, the Id t — J models 

We here reinvestigate the zero temperature transport 
through a LL ring described by the TLM. We do not 
rely on the above mentioned approximation of Ref. and 
present results for arbitrary repulsive interactions, thus 
being closer to situations which can be realized in ex- 
perimental setups. This extends the earlier study and 
clarifies the range of applicability of the central approx- 
imation underlying that work. More precisely, we find 
that the approximation of Ref. [7] leads to qualitatively 
correct results only if the smaller of the two numbers, say 
Pr is either 1, 3 or 5, with the refinement that for p = 3 
or 5 the other number q may not exceed 7. Furthermore, 
the important role of the charge current velocity vj. c as 
a factor accompanying the magnetic flux in all relevant 
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formulas was largely overlooked so far. While, according 
to the approximation made in Ref. [7|, the transmission 
probability exhibits pv F /v s = qv F /v c dips over one pe- 
riod of the magnetic flux, we find that the expression 
vj,c/v s gives a much better estimate on the number of 
dips in almost all cases. It is thus the charge current ve- 
locity and not the velocity v c of the bosonic modes which 
is of relevance for the dip structure of G((j>). We also clar- 
ify how the two expressions for the number of dips and 
their respective partial viabilities are related. 

Transport through LL rings weakly coupled to nonin- 
teracting leads was also studied in Refs. [lj] and [ll| al- 
though with a focus different to ours. 

The rest of this paper is organized as follows. To set 
the stage we introduce the TLM and our way to compute 
the linear conductance in Sec. |TTJ The latter involves 
the one-particle Green function of the TLM which we 
also discuss in this section. In Sec. Mil we describe the 
approximation of Ref. and the corresponding results 
for the flux dependence of the linear conductance. Our 
results for G(4>) are discussed in Sec. IIVI and compared 
to the earlier findings of Ref. |7|. We also comment on the 
relation to the numerical study Ref. H. We conclude with 
a summary in Sec. [Vj 



II. MODEL AND METHODS 

A. The Tomonaga-Luttinger model 

The low-energy physics of Id correlated metals is gov- 
erned by long-range effective forces. Consequently, rele- 
vant wave numbers are from relatively narrow momen- 
tum intervals close to the Fermi momenta ±fcp ^ It 
is then natural to linearize the free dispersion relation 
around the two Fermi momenta. The electrons are sepa- 
rated into two classes, namely into right- and left-moving 
fcrmions, respectively. The Hamiltonian then reads (in 
standard second quantized notation) 

H=H Q + H int , (2) 
H =^ ^av F (k - ak F )c\ a iy c k , a ,<7i 

k a, a 

Hint = ~2L ^a',a'(q) c k+q,a,cr C k' -q.a 1 >' 

k,k' ,q a,a f a,a' 

with a = ± labeling the two branches and a the spin 
direction. The prime at the momentum sums indicates 
that the momenta are restricted to a range 2A around 
the two Fermi points, with the momentum cutoff A > 
0. Since the electrons move in a finite system of length 
L with periodic boundary conditions, the fc-sum is over 
discrete values k n = ^j-n, n e Z. The elements of the 
interaction matrix are conveniently classified according 
to the so-called g-ology convention. 1 — For the TLM the 



following couplings are relevant: 

«£a(«) = 34,11 (?). 5 £«*(9) = 04,.l(?), 

C-a(«) = ^Jl^) , C,-a(9) = 92,±(jl) ■ (3) 

The 34-processes correspond to intra-branch and the gi- 
processes to inter-branch scattering. The || -index in- 
dicates processes involving two electrons of parallel spin 
orientation while _L stands for scattering of electrons with 
opposite spin. 

After normal ordering the operators with respect to the 
ground state and sending the momentum cutoff A to 00 
the model can be solved exactly using bosonization, that 
is the introduction of bosonic creation and annihilation 
operatorsii^ The first step towards bosonization is to 
consider the density operators (q ^ 0) 



Pi 



a,c — C fc,q.cr C fc+g,a,CT ■ 



(4) 



Operators w q a and 6 9)(T obeying bosonic commutation re- 
lations are obtained defining 



27T I Pq,+,<r q > 
\q\L \Pq,-,a q < 



(5) 



These operators as well as the different couplings gi and 
34 may be mapped onto new ones which are no longer 
associated with the two spin directions but with collective 
charge and spin excitations of the system, respectively 
Lv = 2,4) 



9v,c(q) = 9v,\\ +9u,x , 
9v, s {q) = 9v,\\ - 9u,± ■ 



(6) 



(7) 



We define excess particle number operators M a ,ct with 
respect to the ground state (normal ordering). Spin and 
charge operators for a fixed branch index are obtained as 
the linear combinations 



1 



(A^, T ±A^x) , 



charge and spin current operators as (a = c/s) 



(8) 



Ja=N+, a -M-,a (9) 

and the total particle number and spin operators as 

M a =M+, a +N-,a- (10) 

In terms of these newly defined operators and couplings, 
the Hamiltonian splits up into two commuting parts 



H 



TL 



H. 



TL 



H, 



TL 



(11) 
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with 



H, 



TL 



fe>0 



+k 



9iA k ) 
2ir 



( b ljk,a + bl k J- k , a ) 



92,a(k) 



2L 



(12) 



The velocities vj,a an d Ujv.a & r e defined in terms of the 
couplings as 



VN/J,a = VF 



g4 ,q(0) ± g 2 ,q(0) 

2tt 



(13) 



The part of Hj h involving bosonic operators can be di- 
agonalized by a Bogoliubov transformation 1 leading to a 
Hamiltonian of free and uncoupled bosons with disper- 
sion uj a {k) — v a {k)\k\ and velocities 



v a (k) 



VF 1 



9i,a{k) + .g2,a(fc)\ 1/2 



2ltVF 



L | g4,a(fc) ~ 92,a( k ) \ 1/2 ^ 



In the limit fc — + this leads to the charge and spin 
velocities defined in Eq. (fT]). 

In the low-energy limit the details of the momentum 
dependence of the couplings are generically irrelevant as 
long as the g V: \\/±{k) are slowly varying functions close 
to k = (for an exception see Ref . \hm . Without loss of 
generality we thus assume that 



Ju,\\/x(k) = g v ,\\/± ©(fc c - \k\) 



(15) 



with a cutoff k c = 2i:n c / L <C kp of the momentum trans- 
fer. Therefore the velocities v a are independent of k for 
all k < k c . With this choice a closed expression for the 
retarded one-particle Green function Q (x,uS) of the ho- 
mogeneous TLM can be given, which enters the approx- 
imate expression for the conductance through the ring 
used here. To obtain the Green function we follow the 
procedure introduced in Ref. [H (see below). 



It is important to note that the characteristic velocity 
appearing in the flux dependent part of the Hamiltonian 
is the current velocity vj >c . This can easily be seen as 
follows. For a system with an equal number of left and 
right moving electrons (J~ c ) = and the last term in 
Eq. (fl6|) determines the ground state persistent current 
I = —dEo(4>)/d(j), with Eq being the ground state energy. 
It was shown numericall y 1 ^ 18 for lattice models that the 
velocity appearing in the persistent current is the charge 
current velocity vj jC . In Ref. mostly vf was used in- 
stead of i>j, c - Only for Galilean invariant systems both 
velocities are equal^ which indicates that the authors 
of Ref. had such systems in mind without mentioning 
it explicitly. 

We later need to compute the one-particle Green func- 
tion of the isolated LL ring. The term linear in the flux 
of Eq. (fT6)l affects this correlation function. Using Eqs. 
flU) and JSJ the ^-linear term can be written as 



vj. 



-2V2J C 



4> _ vj tC 4> 



2L ' ~^ir L 
x(Ar +iT +jV +4 -AL, T -AL li ) (17) 

and the prefactor of the operator part can be understood 
as an a-dependent correction of the chemical potential 
Sfi a . It thus appears as a phase factor exp (— i S/i a t) in 
the Green function [see Eq. ([20]) ] . 



C. The conductance of the Luttinger liquid ring 

Next, the finite size LL ring is coupled via tunnel bar- 
riers with hopping matrix elements t' to left and right 
leads at positions x = and x = L/2. For simplicity we 
assume equal left and right tunnel barriers. To be specific 
both leads are described as Id tight-binding chains with 
hopping matrix element t, but we expect our results to 
be independent of the precise form of the band structure 
of the reservoirs. 

For noninteracting electrons the transmission probabil- 
ity per spin direction through the ring can be expressed 
in terms of the spin-independent retarded one-particle 
Green function Q r (x,uj) of the isolated ring at position 
and L/2. It is given by£ (for a simple derivation of this 
formula see Sec. Ill of Ref. \2£ 



B. The Luttinger liquid ring with magnetic flux 

We now consider a ring of finite length L described by 
the TLM that is pierced by a magnetic flux <fi. The term 
Hj = irvj tC Jc /(2L) in the Hamiltonian Eq. (|T2|) must 
then be replaced by 



16.17.18 



H 



J 



2L \ Jc V 7T 



J?-2V2J C 



2L 



(16) 



4t 2 sin 2 k 



T( W ) = 



Lt' 2 g H (L/2,to) 



[u - Lt' 2 g R (0, u) + re tk } 2 - Lt' 2 G R {L/2, w) 



(18) 

The zero temperature conductance, on which we focus, 
follows by setting uj — 0, that is the energy equal to the 
chemical potential, and multiplying T by 2e 2 /h. In Refs. 
and d it was argued that this expression can also be 
used in the presence of interactions in the ring provided 
t' is sufficiently small and the Kondo effect does not play 
a role. The latter holds in parameter regimes with an av- 
erage number of electrons on the ring which is even. That 
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the general features of the conductance through an inter- 
acting ring which is not in the Kondo regime arc indeed 
captured to some extent by the approximate Eq. (fl8|) 
was shown in Rcf . using a functional renormalization 
group approach which can be applied for sufficiently weak 
interactions but is nonperturbative in the coupling to the 
leads. Thus we will here also rely on Eq. (fT51) . 



with 



E 



Im+n " -j " J m{l-j)+n ' 



3=0 
I 

,(m) V , 

U lm+n " ■• 



(m) _ \ ^ (l/2m) J ,(m-l) 

| m(l—j)+r. 



3=0 



3- 



(24) 



D. The Green function of the 54-model 

In a first step we study the so-called (74-model in which 
all inter-branch scattering processes are set to zero. We 
decompose the retarded Green function as 

ig R (x,t) = e(t)J2[i^(x,t) + ig<(x,t)] , (19) 

a=± 

with the greater Q> and lesser Q< Green functions of 
right (a = +) and left (a = — ) moving electrons. Using 
bosonization of the field operator—^ for the box-shaped 
potential Eq. (fT5)) . these Green functions are given byi£ 

ig>'<{x,t) = ie i ° t *^ t / L g>/<> (x > t)e F ° /< ^ (20) 
where 



F>/<(x,t) = £i 



_j -^r^aa; — u s t) -^n(aa:— VFt) 



(21) 



The noninteracting Green functions read 



iG>> u (x,t) 



_ iakp-x i ^-(ax—VFt) 

L 



x 
1 

r 



■y _ i -jj-(ax—v F t+iO) 

iakf-x 



i ^-(ax—Vft—iO) 



(22) 



The second exponential factor in the equation for Q> ,0 
appears due to the choice of ±kp — zL2TTnp/L corre- 
sponding to the last occupied level. For L — ► oo it be- 
comes unity, but is relevant in our case as we study rings 
of finite length. As discussed in Ref.llathe Fourier trans- 
form of the greater Green function can be computed it- 
eratively as 



00 ra ra — l 

m=0 i=0 j=0 
(n c ) i(nc)iK) „ia^ 



uj + at(pvj tC /L — 2^ [(m — / — i + 1)i>f + + jVg\ + i0 



(23) 



and initial values 



1] = Y 

a / j 

3=0 

1) _ (1/2)' 



(-Vi) 



Similarly, the lesser Green function reads 



(25) 



L 



00 m m — I 

EEE 

m=0 i=0 j=0 



uj + a<f>Vj iC /L + [( m — I — j)v F + h 



i0 



(26) 



It is easy to see that for \k ± &f| < 9c the one-particle 
spectral function as a function of energy obtained by tak- 
ing the imaginary part of the momentum Fourier trans- 
form of Eqs. (j23|) and (|26|) has support only between 



±v s \k ± ftp] and ±u c |fc ± kp\.— In the thermodynamic 
limit the spectral weight shows a square-root singularity 
at the two edges. The corresponding Green functions (for 
|x|,«j?|t| 3> 1/Qc) hi the (x,t)-plane are given by 



iG>(x,t) 



l 

T 



i^-(a.x — v c t-\-iO) 



-1/2 r 



^ ^i^-(ax — v s t+iQ) 



-1/2 



iG<(x,t) 



L 



2 g — i^jj-(ocx— v c t— i0) 



-1/2 



-1/2 



(27) 



Power-laws with anomalous (interaction dependent) ex- 
ponents only appear if also the gi terms are kept. 1 Below 
we return to Eq. (|27|) when discussing the results of Rcf. 
7. As the anomalous contributions to the propagator 
were neglected in this work, it is the (74-model which was 
effectively studied. 

From Eqs. (fT3")) and (fT4")) it follows that in the ^4-model 
the charge current velocity vj >c and the charge velocity 
v c are equal. 



E. Green function of the full model 

Also for the full model with intra- and inter-branch 
scattering processes a closed iterative expression for the 
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Green function can be given if a box-shaped two-particle 
interaction is assumed. Now, the anomalous dimensions 



la 



s s and 7 C = s c appear. The variables s a are de- 



fined by s a {q) = sinhOg.a at q — and the angle Q q , a 
characterizes via 

(3 q ,a = b q ^ a cosh 6 9:Q - bL q<a sinh e 9jQ (28) 

the canonical transformation to new annihilation (cre- 
ation) operators /3 q ]l by means of which the bosonized 
TL-Hamiltonian is diagonalized. The angle Q q , a depends 
on the couplings according to the formula 



tanh(29 g , a ) 



2nvp 

9i,\\(q)±g4,±(q)^- ] 



2nvp 



(29) 



where the "+" -signs are relevant for the charge angle Q q . c 
and the "—"signs for the spin angle O g , s - The function 
F> f< (x,t) in Eq. (20]) now reads 



c i i 

^ >7< m = E~[ 



— e ^ 



e ±i^-n{ax-v F t) g 2 _ g 2 



n=l 

2 1,2 

_|_ ^gT'T™^^') 4- s c ±i^n(ax-v s t) j-gg-j 

.2 

i *s ^i^-n(ax+v s t) 

+ 2 

Because of the assumed box-like shape of the two-particle 
interaction (in momentum space) the s a (q) are indepen- 
dent of q = q n = 2Tm/L for all n < n c . 

We again introduce recursively defined coefficients. 
The at 1 ' are given as in Eq. and coefficients bm°a 
as well as Cm^l are computed according to the formulas 



dm) 



a ^ y 



3=0 



([l+ 7a ]/[2m]) J (m _!) 

j\ m(l-j)+i,a ' 



I 



c (m) = > 

3=0 



(7a/[2m]) J (TO _!) 



(31) 



-'m(i— j)+z,a 



with initial conditions 

_ ([l+7a]/2) T 



.(1) 



to! 

(7a/2) r 



(32) 



In analogy to Eq. (|23|) one can bring the greater and 
lesser Green functions into the form 



9>(x,u>) = A 



- 4 -7c -7s 



ia(k F + 2jj-)a 

L 



oo m I m—l p 



E E E E E ^.fife^c^ft^ 



m=0 i=0 j=0 p=0 g=0 



m—l—pl—j,cj,c "p—q,s^q,. 



e ia^-(m-2j-2q)x 



uj + ac])vj c /L — [[(m + 1 — I — p)vf — lv c — pv s ]) + «0 

(33) 



and 



a iakpx 



g<{x,w)=A-T 

oo m i m—l p 

EEEEEeUfe^et^ 

m=0 i=0 3=0 p=0 g=0 

g — ict2f(m— 2j— 2g)x 

+ a4>v,;_ c /L + ^ ([(to - / - p)v F - lv c - pv s ]) + i0 

(34) 



III. THE APPROXIMATION OF JAGLA AND 
BALSEIRO 

As mentioned in Sec. HH in their work Jagla and Bal- 
seiro (JB) 7 effectively studied the ^4-model as they ne- 
glected the anomalous contributions to the propagator. 
The starting point is the retarded Green function given 
by the sum of the approximate expressions (|27[) (valid 
for |ar| , 1 1 1 ^> l/q c ), where in addition vj tC is replaced 
by vp. These two velocities become equal in Galilean in- 
variant systems with g^ c (Q) = <72.c(0) (and trivially in the 
noninteracting case). This condition cannot be achieved 
within the (74-model for any reasonable choice of 11 and 
<74,_i_. It is thus to some extend inconsistent to replace 
the charge current velocity by the noninteracting Fermi 
velocity. Compared to Eq. (|2T|) . JB also seem to have 
neglected the factor exp {i ^(ai — vpt)} appearing in 
G>(x,t) [see Eq. (6) of Ref. 7] which has to be included 
due to the finitcness of the system. 

To analytically perform the Fourier transformation of 
the approximate Green function from time to frequency 
JB only consider cases where v c /v s — q/p, with q and p 
being "small" odd numbers, which corresponds to rather 
specific choices of two-particle couplings. They assume 
that the Fourier transform is essentially determined by 
the behavior of the "dominant poles". For x = L/2 
poles of the form l/(t — t*) appear at t* — qL/(2v c ) — 
pL/(2v s ), that is, times at which charge and spin exci- 
tations starting (at the same time) from the left contact 
and traveling along the ring with v c and v s get together 
at the right contact. One expects this approximation to 
work best for small numbers p and q because otherwise 
such poles are rare among all poles of G> (x,t). The func- 
tion 



e(t) 



E 

a=±l 



e i^j-(ax-v F t) 



e e 



■ (ax— pv c t-\-i0) 



+ 



1 



- (ax— pv c t— iO) 



(35) 



has for x — and x = |- the same "dominant poles" 
i* and residua as the Green function resulting from sum- 
ming up the two terms of Eq. |27|) . The time between two 
such poles is given by T\ — ^f-T$ = ^-Tq. The time 
Tq = vp I L is the <j> = period (in t) of the Green function 
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in the noninteracting case, while T\ is the corresponding 
period of the approximate expression (|35[) including the 
interaction. Fourier transforming Eq. (|35|) leads to 



J-1^V C V S a=±1 



iakpx 



E 



uj + acj)vj, c /L - *z(m + l)jr- + iO 



-iot-r-mx 



uj + axf>Vj tC /L + ^m±r + iO 



(36) 



Within this "dominant pole" approximation and for 
v J,c — ► vf, the shape of the transmission probability Eq. 
([T5| depends almost exclusively on Ti/T . Replacing T 
by Ti, mainly rescales the </>-axis by a factor Tq/Tx. In the 
noninteracting case T(</>) is a periodic function of period 
2tt which is symmetric with respect to <j) — -k. Replacing 
<p by 4> + 27rm To / T\ so that the new magnetic flux is still 
within [0, 2f\ the transmission probability remains invari- 
ant as one has effectively shifted the summation variable 
ui — ► m + n in Eq. (f36|) . Thus, the sole effect of the 
different interaction dependent velocities is to reduce the 
periodicity of the magnetic flux dependence of the con- 
ductance by a factor Tq/T\. Since in the noninteracting 
case the number of dips of G(</>) = e 2 T(</>) /h in the inter- 
val of periodicity equals 1 (see Fig.[TJ), it becomes Ti/Tq 
in the presence of interactions according to the "domi- 
nant pole" approximation. This effect is shown in Figs. 
2 and 3 of Ref. 0. 

To some extend the "dominant pole" approximation is 
motivated by the simple physical picture that the spin 
(traveling with v s ) and the charge excitations (traveling 
with v c ) have to "recombine" at x = L/2 (and x = 0) for 
an electron to pass the LL ring. Indeed, the poles kept 
in the present approximations correspond to such times. 
We next show that the usefulness of this picture is quite 
limited and that the number of dips of the conductance 
resulting from an exact evaluation of the Green function 
differs considerably from the above result even in the 54- 
model. 



dependent on the length of the ring L, provided L is 
not too small. Calculations in this paper are performed 
for L = 256vf which turns out to be sufficiently large 
to observe structures as described in Ref. 0. To bring 
out well discernible dips in the transmission probability 
and thus results comparable to those obtained by Jagla 
and Balsciro, who neither give their t' nor the size of the 
energy window averaged over, we choose t' 2 = 0.005 and 
take the average over the 10- interval [— w max , Sa] = 
[—0.05, 0.05] if not mentioned otherwise. The width of 
the box potential is determined by n c , see Eq. (|15p. which 
we choose to be n c = 5 but since we are dealing only 
with low-energy properties of the system, all results are 
practically independent of the width of the box potential, 
as we have also checked numerically. 




FIG. 1: Transmission probability for the noninteracting case 
in which all characteristic velocities become equal to vf- 

Figure [T] shows the transmission probability as a func- 
tion of the magnetic flux in the noninteracting case, i. e., 
in the case where all characteristic velocities are equal to 
vp- It displays a "dip" at (f> = it which occurs indepen- 
dently of the chosen coupling to the leads and the aver- 
aging frequency interval [— w max , Sn]' Varying these 
parameters either renormalizes the curve as a whole or 
makes the dip sharper or wider. 



A. Results for the 34-model 



IV. RESULTS 

We now discuss the conductance (transmission) 
through the ring described by the TLM and coupled to 
leads at positions x = and x = L/2. We use the ap- 
proximate expression (|18p relating the retarded Green 
function and the transmission probability. 

The transmission probability is averaged over a small 
energy window around u> = 0. This has also been done in 
Refs. @ and [H Hallberg et al& argue that such an aver- 
aging over a finite energy window accounts for "possible 
(gate and bias) voltage fluctuations and temperature ef- 
fects unavoidably present in an experimental system" . 

Our energy unit is fixed by setting the hopping in the 
leads to r = 1. Transmission curves are only weakly 



We next compute the transmission probability Eq. (|18p 
as a function of <f> using the exact Green function [see Eqs. 
(HHJ), (H2J), and <HH])] of the (74-model with a box-shaped 
two-particle potential. 

We first consider cases where the interaction param- 
eters are chosen such that v c /v s = q/p with odd inte- 
gers q and p, as well as vj^jvp — 1, i.e. cases where 
the point made about which velocity to use when taking 
into account the magnetic flux is irrelevant. This cor- 
responds to the situation mainly studied by Jagla and 
Balseiro. We again emphasize that vj >c = v c for the 
<74-model and the charge velocity and the charge cur- 
rent velocity can be identified throughout this subsec- 
tion. Our calculations show that under the above condi- 
tions the "dominant pole" approximation remains a good 
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guide only for Ti/Tq ;$ 7 but that it loses its validity for 
Tx/Tq 2; 9. In particular, for sufficiently small Ti/T 
the number of prominent dips within [0, 2-7r] is given by 
Ti/Tq — qvF/v c — pvf/v s . This is shown in Fig.[5]where 
curves for different values of Ti/Tq are presented. While 
for all cases where T x /T = 5 [Fig. (a)] or T x /T = 7 
[Fig. [2] (b)] the number of prominent dips does indeed 
equal 5 or 7, respectively, it is 9 for Ti/Tq — 9 only if 
p = 1, i.e. v c /v s = 9 [solid line in Fig.[2](c)]. The exam- 
ple of the curve with v c /v s = 9/5 [dashed line in Fig. [5] 
(c)] shows that in this case the numbers p = 5 and q = 9 
are apparently not small enough but that the transmis- 
sion probability rather resembles that of the case where 
Vc/vs = 5/3 [dashed line in Fig. [2] (a)]. Further down 
we argue that the latter is not accidental. Our analy- 
sis specifies how small the odd integers q and p have to 
be for the "dominant pole" approximation to be applica- 
ble. Specifically, we find that this approximation leads to 
qualitatively correct results only if the smaller of the two 
numbers, say p, is either 1, 3 or 5, with the refinement 
that for p = 3 or 5 the other number q may not exceed 
7. 



We now proceed and consider situations in which the 
coupling constants are still chosen such that v c /v s = q/p, 
but v,j^ c /vf = v c /vf > 1 (for the 54-model with repul- 
sive two-particle interaction). The factor Vj tC in front 
of the magnetic flux in Eq. (|36[) . leads to a decrease of 
the periodicity of T(</>) by a factor vf/vj^ c and thus to 
an increase of the number of dips by a factor vj, c /vf 
with respect to Ti/Tq, the result obtained for vf = vj >c . 
For sufficiently small q and p, that is if the "dominant 
pole" approximation is applicable, we thus expect to find 
(vj,c/v F )(Ti/T ) dips. This is illustrated in Fig. [3] for 
vj jC — v c — 3Ui?/2 = 3v s , that is p = 1, q = 3. 
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FIG. 2: (Color online) Transmission probability for the g^- 
model with vj, c /v a = v c /v s = q/p, vj, c /v F = 1 and Ti/Tq = 
qvF /v c = pvf /v a = 5, 7, 9 from top to bottom. The number 
of prominent dips is always equal to T\ /To with the exception 
of the case where v c /v 3 = 9/5 (see the text). The coupling 
to the leads is t' 2 = 0.001 in these pictures to minimize the 
dependence of the curves on the energy interval over which 
the average is taken. 



FIG. 3: Transmission probability for the (74-model with 
v c /v s = 3/1, vj. c /vf = v c /v F = 3/2 > 1 and 
(«j, c /«f)(Ti/T ) = 3 dips, while Ti/T = 2. 

If the "dominant pole" approximation were applicable 
for all p, q, one would thus have (vj^ c /vf)(Ti/To) dips in 
T(4>) in the interval <f> 6 [0, 27r]. Given any arbitrary com- 
bination of the (relevant) velocities v s , v c , i>j, c , a natural 
strategy for guessing the number of dips would be to look 
for "small" odd numbers p and q for which pv c « qv s and 
to compute (vj, c /vf)(Ti/T ) from these. The integers p 
and q would probably be looked for so that pv c s» qv s 
holds as accurately as possible for, at the same time, p 
and q as small as possible. From the dashed line in Fig. [2] 
(c) it is clear that already the numbers p = 5 and q = 9 
are not small enough for the "dominant pole" approxi- 
mation to be reliable. In such cases, one has to resort 
to smaller odd numbers p and q — for which pv c « qv s 
holds less well — to predict the number of dips. As al- 
ready shown in Fig. [2] (c), for v c = vp and v c /v s = 9/5, 
instead of p — 5 and q — 9 one has to take p = 3 and 
q = 5 to obtain the correct number of dips. As we shall 
further show, it is mostly necessary to choose p and q 
such that for the smaller number p — 1 even if the ap- 
proximate identity pv c w qv s does then hold only to a 
very low degree of accuracy. This leads to our central 
result 



vfT 



p=i 



V.J,cVF 

v F v s 



Vj,c 



(37) 



for the number of dips. 

Consider, as an example for the validity of Eq. (|37p . 
the case where v c /vf — H^s/(3t | F) = 11/6 (p = 3, 
q = 11). We then have Ti/T = pv F /v s = 3-2 = 6. 
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FIG. 4: Transmission probability for the <?4-model with the 
ratio of charge and spin velocities as indicated. In all cases, 
the number of prominent dips is given by one of the integers 
closest to vj, c /v s (= Vc/vs for the <?4-model). 



Taking into account that Vj tC /vF = v c /vf = 11/6, follow- 
ing the "dominant pole" approximation one would arrive 
at the conclusion that the number of dips should equal 
vj t cT\/ '(vfTq) = 11/6-6 = 11. In fact, the numbers p = 3 
and q = 11 are already too large for this approximation 
to be valuable. Smaller numbers, for which the equation 
pv c = qv s is still "approximately" fulfilled, are p = 1 
and <? = 3. This gives Ti/T | p=1 = pv F /v s \ p=1 = 2 and 
thus (vj, c /v F )(Ti/T )\ p=1 = 11/3. The number of dips 
should, of course, be an integer number, and the result 
11/3 lies between 3 and 4 so that one can expect one of 
these two to be the right answer. Indeed, as shown in 
Fig-Si (a), the corresponding curve exhibits 3 dips. This 
kind of reasoning is applicable whenever the approximate 
identity pv c « qv s does not hold to some sufficiently high 
degree of accuracy with p = 3, 5 and q — 5,7 [as it was 



the case for the example of the dashed line in Fig. [2] (c)]. 
We find that Eq. ([3T| for estimating the number of dips 
of T((j>) can be used in almost all cases. Further examples 
of its validity are shown in Figs. 2] (b), (c) and (d). For 
Fig. [4] (b) vj^c/vg — 3.8 and the number of dips is four, 
while for Fig. [4] (c) v c /v s — 25/2 so that the numbers p 
and q for which pv c — qv s are no longer both odd. Equa- 
tion (|37p . however, applies nevertheless and, in this case, 
leads to the result vj tC /v s = 12.5 again giving a correct 
estimate of the number of dips. 

Jagla and Balseiro argue that, according to the as- 
sumptions ("recombination" of spin and charge) under- 
lying their "dominant pole" approximation, the transmis- 
sion probability would have to be "very small" whenever 
the couplings are chosen such that {v c /v s ) ±1 — 2n with 
integer n. In these cases no poles of the form l/(t— t*) 
appear in Q(x = L/2,t) and charge and spin excitation 
cannot "recombine" at finite times at x — L/2. 7 However, 
the expectation that the transmission probability is small 
under these conditions is not confirmed by our calcula- 
tions. Fig. [4] (d) shows the example of a curve for which 
v c /v s = 4/1 and in this case (and in all the others we 
have checked) the transmission probability is not partic- 
ularly small. This observation provides further evidence 
that the idea that charge and spin excitations have to 
"recombine" at x = L/2 for an electron to pass through 
the ring is of very limited usefulness in the present con- 
text. 



B. Results for the full model 

All the results obtained in the model with only <? 4 - 
couplings essentially carry over to the full model (with 
additional inter-branch interaction g^), the main differ- 
ence being that now u,/. c 7^ v c . The charge current veloc- 
ity vj :C Eq. (fl~3|) becomes smaller than the charge velocity 
v c due to the presence of nonvanishing <?2-couplings and 
it is even reduced to the Fermi velocity if 52,||/±(0) — 
54.||/-l(0)- An example of such a case is shown in Fig. 
[5] (a) and the number of dips agrees with the prediction 
made from Eq. (JSZD vj,c/v s = 1/0.21 « 5. 

In case all four couplings g u ,\\/± are chosen to be equal, 
also the spin velocity becomes equal to the Fermi velocity, 
one then has v s = vj iC = vf- The transmission probabil- 
ity for such a case is presented in Fig. [5] (b). It shows a 
dip structure similar to that of the noninteracting case. 
This, however, does not come as a surprise in view of 
the foregoing considerations but is in accordance with 
Eq. (33) vj. c /v s = 1. 

Finally, the transmission probability for a more general 
case is presented in Fig.[6l Here, T(<f) exhibits 10 dips 
and this does indeed represent the ratio of the charge 
current and the spin velocities following from the chosen 
coupling constants, namely vj_ c /v s « 1.95/0.195 = 10. 

Our general finding that in most cases vj^ c /v s gives an 
estimate of the number of dips, may be compared to the 
result of Hallberg et al. s . They studied the Id t-J model. 
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FIG. 5: Transmission probability for the full model and dif- 
ferent sets of coupling parameters. The individual couplings 
are chosen such that the (relevant) velocities are vj <c /vf = 1, 
v c /vf ~ 1.79, v s /vf ~ 0.21 in (a) and Vj >c /vf = v s /vf = 1, 
v c /v F » 2.01 in (b). 




FIG. 6: Transmission probability for couplings such that the 
(relevant) velocities are vj, c /vf ~ 1.95, v s /vf ~ 0.195, and 
Vo/vf » 2.32 leading to vjj/v s w 1.95/0.195 = 10. 



The Green function of a few lattice site (up to 8 sites) iso- 
lated t-J ring was computed using exact diagonalization. 
For J = this model corresponds to the U — > oo - limit 
of a one-dimensional Hubbard-Hamiltonian for the ring. 
Whenever the intermediate state contains N electrons, 
the transmission probability shows N equidistant dips. 
The authors argue that in the limit J — the ratio of 
the spin and charge velocities becomes v c /v s = N and 
this way try to relate their findings to the ones of Ref. lZ|. 
It remains, however, unclear how this would explain the 
observed dip structure in view of the results presented 
here which unambiguously suggest that vj tC and not v c 
is the velocity of importance for the number of dips. In 
particular, in the U — > oo - Hubbard model, one finds 
vn^ c = 2v c (see Ref. [l]) so that vj tC — l/2v c ^ v c [see Eq. 



©] 



V. SUMMARY 

In the present paper we have studied the effect of the 
different velocities characterizing the low-energy physics 
of a LL on transport through a one-dimensional, metallic 
ring of correlated electrons. A magnetic flux piercing the 
ring was included and the dependence of the transmission 
probability (linear conductance) on the flux was studied. 
It was obtained using a relation [Eq. Q18p] between the 
one-particle Green function of the isolated ring and the 
transmission probability. This becomes exact in the non- 
interacting limit and serves as an approximation in the 
presence of two-particle interactions. Results were av- 
eraged over a small energy window around the chemical 
potential just as in Refs. andS It proved possible to 
improve on an approximation suggested in Ref. and to 
exactly calculate the Green functions of the isolated ring 
for the case of a potential with box-like shape in momen- 
tum space. Since only the low-energy properties of the 
system are relevant, this restriction to a specially shaped 
potential should not be regarded as severe, in particular 
because the width of the potential plays practically no 
role for the results. 

Characteristic dips appear in the transmission proba- 
bility as a function of the magnetic flux whose number 
depends on the different velocities in play. According to 
the "dominant pole" approximation of Ref.[7|, the num- 
ber of dips for (j) G [0, 27r] is given by pvf/v s = qvpjv c 
(or pvj^ c /v s — qvj^ c /v c if, as for systems which are not 
Galilei invariant, Vj tC ^ vp) if small odd integers p and 
q can be found such that pv c — qv s . For generic repul- 
sive two-particle interactions v s < v c and p will be the 
smaller of the two integers. We have shown that p and 
q must indeed be very small for this prediction for the 
number of dips to hold. For more generic choices of the 
relevant velocities v c , V s , and vj tC , and thus coupling con- 
stants, the number of prominent dips in T(<f>) over one 
interval of periodicity can be estimated to be an integer 
close to Vj tC /v s . The wide applicability of this relation 
was demonstrated. 

That the charge current velocity vj tC is a relevant pa- 
rameter determining the flux dependence of the linear 
(charge) conductance through the ring should not come 
as a surprise. The appearance of the velocity of the collec- 
tive spin excitations v s can be understood from the struc- 
ture Eqs. (23J), (EH), ©, and (J3U) of the one-particle 
Green function. Generically the sums appearing in these 
expressions are dominated by the terms in which the in- 
teger in front of v c is zero while the integer in front of v s 
is greater than zero (but small). 

A simple picture related to the "dominant pole" ap- 
proximation suggests that electrons can pass the ring 
only if spin and charge excitations traveling with different 
velocities (spin-charge separation) can recombine at the 
right contact. Our results show that this idea of charge 
and spin "recombination" is of limited usefulness for the 
understanding of transport through LL rings. 
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